Nano-pattern induced ferromagnetism in strongly correlated electrons 
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Band ferromagnetism in strongly correlated electron systems is one of the most challenging issue in today's 
condensed-matter physics. In this theoretical work, we study the competition between kinetic term, Coulomb 
repulsion, and on-site correlated disorder for various lattice geometries. Unconventional and complex ferromag- 
netic phase diagrams are obtained: wide region of stability, cascade of transitions, re-entrance, high sensitivity 
to the carrier concentration and strongly inhomogeneous ground states for relatively weak on-site potential. The 
direct and systematic comparison with Exact Diagonalization shows that the Unrestricted Hartree-Fock method 
is unexpectedly accurate for such systems, which allows large size cluster calculations. A match of the order 
of 99.9% for weak and intermediate couplings is found, slightly reduced to about 95% in the large repulsion 
regime. Nano-patterned lattices appear to be particularly promising candidates that could, with the tremendous 
progress in growing and self-organized techniques, be synthesized in a near future. 
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The possibility of ferromagnetic ground state (GS) in 
strongly correlated itinerant systems is a highly non triv- 
ial question and is still one of the most debated major is- 
sue. The minimal Hamiltonian, introduced fifty years ago by 
Gutzwiller, Hubbard and Kanamori in order to tackle such a 
subtle physics, is known as the Hubbard model 1 . More pre- 
cisely, one of the main issue was the better understanding of 
the appearance of ferromagnetism in transition metals. The 
Hubbard model was introduced fifty years ago by Gutzwiller, 
Hubbard and Kanamori in order to get a better understanding 
of the appearance of ferromagnetism in transition metals^. A 
first rigorous example of ferromagnetism was provided by Na- 
gaoka and Thouless^ a couple of years later; a certain class of 
Hubbard models exhibit saturated ferromagnetic GS for infi- 
nite Coulomb repulsion at half filling and under the condition 
of having exactly one single hole in the system. The question 
of whether Nagaoka ferromagnetism survives at finite repul- 
sion and at larger hole density has not been answered so far. 
More than twenty years later, Lieb established a theorem 4 , stat- 
ing that a system of half-filled bipartite lattices, with different 
numbers of sublattice sites, possesses a unique ferromagnetic 
GS (in fact ferrimagnetic). In the same period, another im- 
portant rigorous result was obtained by Mielke 5 and Tasaki 6 . 
They formulated a theorem stating that systems with nearly 
fiat and partially filled band have a global stable ferromag- 
netic GS. For example, a very beautiful striking example of 
this kind of ferromagnetism is achieved on the Kagome lattice 
(3 atoms per unit cell). It is worth noticing that Lieb ferrimag- 
netism appears to be a particular case of Mielke-Tasaki's mag- 
netism at half filling. One remarkable feature of this theorem 
is that even an infinitesimal coulombian repulsion U leads to 
a fully spin polarized GS while the magnetism is induced by 
the electronic correlations. More recently, several numerical 
studies have shown that ID lattices based on connected trian- 
gles are favourable for ferromagnetis m 7 ' 8 ' 10 . This emphasizes 
the important role of the frustration in the statibilization of fer- 
romagnetism. In the same spirit, the possibility of ferromag- 
netism in quantum dot super-lattices 11 and organic polymers 
consisting of chains of five-membered rings 12 have been in- 
vestigated leading to possible experimental realizations as the 
polyaminotriazole 13 or polymethylaminotriazole 14 . 
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FIG. 1: (Color online). The three geometries of the nano-patterned 
lattices considered in this work: (a) square, (b) triangular and (c) 
graphene. The defect plaquette (shadded polygons) sublattice is de- 
fined by the elementary vectors (arrows). 



Another exotic type of ferromagnetism has been recently 
reported in oxide thin films, e.g. HfOa 15 and CaO 16 , also 
known as d° ferromagnetism. In this case, the presence of non 
magnetic defects (either vacancies, cationic and/or anionic 
impurities) induces finite magnetic moment in their vicinity 
and eventually long range ferromagnetic order. These unex- 
pected findings have been followed by numerous debates^—. 
A simple interpretation of this d° ferromagnetism has been 
proposed within a single orbital Hubbard model with short 
range correlated disorder on the oxygen orbitals by using an 
unrestricted Hartree-Fock method 20 . This study also led to 
interesting predictions; the compound A 1 _ K B ;1 ,02, A being 
either Zr, Ti, Hf, etc, doped by K or Na, could provide high 
critical temperature ferromagnetic systems, as confirmed later 
on by ab-initio based calculations 21 . Finally, this scenario 
has also explained the possibility of ferromagnetic phases 
in graphene and irradiated graphite 22 ' 23 . The recent devel- 
opments in nanotechnology growth and self-organized tech- 
niques are opening new paths towards the realization of quasi 
ID and 2D super-lattices, nano-cluster arrays 24 or decorated 
lattices 25 . Hence, ferromagnetism in strongly correlated elec- 
tron systems is still an uncharted terrain worth to be explored. 

In the above mentioned numerical studies, the crucial in- 
gredient is the correlated nature of the disorder. For exam- 
ple, in K doped ZrC»2, the substitution of the cations Zr + by 
K + provides three holes per impurity and induces short range 
correlated onsite disorder on the neighboring oxygen orbitals. 
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Thus, the main motivation of this study is to analyse the ef- 
fects of short-range correlated disorder in two dimensional in- 
teracting electron systems. More precisely, we investigate the 
possibility of ferromagnetism on various nano-patterned lat- 
tices. In order to draw general conclusions, we have consid- 
ered three different geometries: square, triangular and hon- 
eycomb (graphene) lattices, as depicted in Fig. Q] Note that 
the origin of such defect plaquettes in real material could be 
the substitution, presence of add-atoms or vacancies, or other 
structural defects. 

The minimal model for that purpose is the one-orbital Hub- 
bard Hamiltonian : 

<i,j> i z£def..cr 

where c\ — (c|^., c]j_) is a spinor of fermion (electron) cre- 
ation operators at site i and = c\ a Ci, a the on-site den- 
sity operator. Note that the hopping term t and the on-site 
Coulomb repulsion U have the same magnitude on all bonds 
and sites respectively; the on-site spin independent scattering 
term V is defined only on sites belonging to certain defect pla- 
quettes of the nano-patterned lattice as depicted in Fig. Q] In 
what follows, U and V will be expressed in units of t. De- 
fect plaquette sublattices on the three geometries considered 
in this work are also displayed. These super-lattices are de- 
fined by unit-cell vectors (ti, t 2 ), both expressed in the units 
of the underlying lattice vectors. For example, in Fig.[T}b (tri- 
angular lattice), these vectors are respectively ti = (3,1) and 
ta = (0,3). 

To handle such kind of models, various numerical meth- 
ods have been employed during years, all with advantages and 
limitations. On one side, exact diagonalization (ED) allows an 
exact treatment but is restricted to small clusters only due to 
an exponential growth of the Hilbert space. On an other side, 
within Quantum Monte-Carlo (QMC) based methods, larger 
sizes are reachable, but the sign problem can occur. Finally, 
density matrix renormalization group (DMRG) can circum- 
vent these limitations, but remains essentially suitable for ID 
systems. 

From the mean field approaches, it is well known that the 
standard (restricted) Hartree Fock approximation, which leads 
to the S toner criterion (Upp > 1, with pp the single particle 
density of states at the Fermi level), overestimates the ten- 
dency toward ferromagnetism in most systems^. It is worth 
noting that in Nagaoka, Lieb and Mielke-Tasaki ferromag- 
netism, this criterion is always fulfilled since pp = oo. How- 
ever, the unrestricted Hartree Fock approach (UHF) has been 
successfully employed in the description of some exact GS 
properties in the presence of disorder. In this approximation, 
the interaction term is decoupled as 

"<i,t n U (™i,t) n *4 + ™i,t( n a) ~ ( n «,tX n u) ( 2 ) 

and the set (rii M ), for i € [1, N] with N the number of sites, 
is calculated self-consistently. It is worth noticing that such an 
approach is, due to the self consistency, non-linear and non- 
perturbative in essence; for a 2D square lattice at half filling 

for example, one gets a gap of the form Au ~ e~ 2?I V t / c/ i 



as obtained in 9 . As a first example, a fairly good agreement 
between this method and ED has been obtained in the case 
of persistent currents of disordered mesoscopic rings 2 -. The 
comparison with QMC also supported the reliability of the 
UHF for the study of disordered superconducting systems 27 . 
In these works, the localized and/or inhomogeneous nature of 
the one particle wave functions is at the origin of the fact that 
the Slater determinant becomes a good estimate of the exact 
many -body GS. It is worth noticing that the UHF usually pro- 
vides a reasonable order of magnitude of the physical observ- 
ables and a rather good qualitative description of phase dia- 
grams. Surprisingly enough, in spite of the expected richness 
of the physics, the validity of the UHF approach has not been 
provided so far for multi-band/orbital systems. We also ex- 
pect in the case of super-lattices (several atoms per unit cells), 
giving rise to complex multi-band structure, the presence of 
relatively flat band regions which should favor the appearance 
of ferromagnetism. Our goal is now to investigate such a pos- 
sibility in nano-patterned lattices. Simultaneously, we will es- 
tablish the accuracy of the UHF approach by direct systematic 
comparison with ED calculations up to the largest accessible 
system (e.g. 8x8 sites). In addition, it should be stressed that 
the agreement between ED and UHF is not only achieved for 
weak and intermediate Hubbard repulsions, but also for large 
U. This will be shown in what follows. 
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FIG. 2: (Color online). Spin resolved density of states on a 72 x 
72 site cluster for respectively (U, V) = (0, 0) (a) and (U, V) = 
(1.5, 1) (b). The carrier concentration is n = 1.5625. The vectors 
characterizing the pattern are ti = (2,2) and t2 = ( — 2,2) (see 

Fig.ru. 

Let us first start with the effects of nano-patterns on the 
one-particle density of states (DOS). For that purpose, we per- 
formed the UFH calculation on a 72 x 72 triangular site cluster, 
at a carrier concentration close to the Van Hove singularity of 
the clean case, namely n = 3/2 + e, with e = 0.0625 (see 
Fig|2). The (U, V) parameters are set to (1.5, 1). As can be 
seen, for finite U and V, the ground state is ferromagnetic. 
Beside this band splitting, one clearly sees in the vicinity at 
the Fermi level, that the Van Hove singularity is now replaced 
by a double-peak structure. This feature survives in the ther- 
modynamic limit. Unfortunately, the DOS is not the appro- 
priate quantity to provide any information on the nature of the 
many body ground state. Furthermore, as pointed out before, 
the standard HF is known to overestimate the tendency to fer- 
romagnetism. Thus, one of the crucial questions is whether or 
not the true GS is really ferromagnetic? In other words, is the 
UFH approach reliable for such systems? As it will be seen 
in what follows, even for such relatively small parameters, the 



3 



GS turns out to be very inhomogeneous with strong fluctua- 
tions in local charge and spin densities, not reflected by the 
DOS. 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 



u u u 

FIG. 3: (Color online). Ground state phase diagrams for the square 
lattice obtained by ED (symbols) and UHF (continuous lines) meth- 
ods. The number of sites N = I x I with I, the carrier (electron/hole) 
number and the total spin S are indicated. In each case, we consid- 
ered 2 defect plaquettes separated by a vector (2,2) in units of their 
respective lattice vectors (see Fig.fT}. 

In order to reveal these drastic effects, we propose to pro- 
ceed with a systematic comparison between ED and UHF cal- 
culations. We start with the ([/, V) phase diagrams for both 
electron/hole doped square and triangular lattices. The results 
are plotted in Fig. [3] and Fig. [4] respectively. In all cases, the 
systems contain two defect plaquettes. Phase diagrams on the 
square lattice for both electrons and holes at filling n = | 
in Fig. 0a)-(b) look identical, as expected for systems pre- 
senting the electron/hole symmetry (electron <H- hole leads to 
V <-> —V). First of all, one immediately sees a remarkable 
agreement between the ED critical lines and the UHF ones. 
We find an extended large spin GS 5 = 2 sector (the largest 
possible spin being 5 = 3 in this case), separated from a sin- 
glet region (5 — 0) without any intermediate triplet phase. 
The critical line of this phase transition is well approximated 
by V c ~ J7/4, hence showing that for a given V, a reason- 
able Hubbard repulsion is enough to stabilize a ferromagnetic 
phase. 
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FIG. 4: (Color online). Ground state phase diagrams for the trian- 
gular lattice obtained by ED (symbols) and UHF (continuous lines) 
methods (see Fig[3]for notations). 

We now consider, for the same carrier density and same 
cluster size, the case of the triangular lattice depicted in 
Fig. SJa)-(b). The phase diagrams appear to be more com- 
plex and richer. In contrast to the square lattice, and as a 
consequence of the particle/hole symmetry breaking, they are 
now completely different. In the case of electron doped, for 



a given U, a cascade from 5 = 2 sector to the singlet one 
through an intermediate 5 = 1 region is obtained. However, 
the hole doped case behaves differently. We indeed observe a 
wide fully polarized 5 = 3 region, and for a fixed value of U, 
three successive phase transitions at V c i, V c ,2 and V Cf 3. These 
critical values correspond respectively to magnetic transitions 
Si — >• 5j of —> 2 —> 3 — > when tuning V from negative 
to positive values. Remark that a fourth additional transition 
between 3 — >• 1 sectors at V > is obtained within the UHF. 
This implies the presence of a 5 = 1 narrow region which 
we interpret as a consequence of the quasi-degeneracy of the 
5 = and 5=1 GS energies, as verified in the ED spectrum. 
However, the agreement between the two methods remains ex- 
cellent in both cases of electron and hole doped systems and 
surpasses our initial expectations. Let us now consider larger 
cluster sizes. Results for the 6x6 site cluster filled with 4 
electrons (n = 1/9) are depicted in Fig. Oc) and Fig. HJc) 
for respectively square and triangular lattice. In contrast to 
panels (a) and (b), the phase diagram plotted in (c) is now 
completely asymmetric. It is not useless to point out that the 
4x4 and the 6x6 systems are inequivalent. For a given U, the 
5 = 1 phase stability extends to very large positive values of 
V, while for V < 0, the triplet-to-singlet transition occurs at 
V c ~ — U /2. The triangular lattice case is even more remark- 
able; the triplet phase stability region is also very extended, 
with a very steep slope for the critical line which could be 
approximated by V c ~ U/4 for repulsive V. But the most 
surprising feature appears in the attractive V region where a 
spectacular re-entrance of the triplet phase for 0.5 < U < 1 
is observed, showing the complexity of the band structure on 
frustrated geometries with the presence of such defect plaque- 
ttes. Once again, even with such a complexity, the UHF ap- 
pears to be in perfect quantitative agreement with ED results, 
hence benchmarking such a mean-field approach for deco- 
rated lattices. Finally, Fig. [3] and Fig. |4] nicely illustrate that 
electron-hole symmetry breaking, complex multi-band struc- 
ture and geometrical frustration are essential ingredients in or- 
der to stabilize and strengthen magnetic phases 7 ^—. 

We now proceed further by focusing our attention on the 
nature of the GS. For that purpose, we analyze the local charge 
and spin densities defined as 

Pi = (ni,f) + (n a ) (3) 

s . z = M - (go) (4) 

for both the triangular and graphene lattices. The results are 
depicted in Fig. [5] and Fig. [6] Prior to this study, we have 
computed within the UHF the GS phase diagrams for each 
system in order to determine interesting sets of (U, V) param- 
eters; chosen to provide a triplet GS in both cases FigJ5}a) and 
Fig|6]and a quintet in case Fig|5jb). Note that for the 8x8 
system (b) we considered only parameters for which the GS 
is fully polarized. Indeed, this spin sector is the only reach- 
able one by ED for such a large system size. Back to the 
left panels of Fig. [5] where a systematic comparison between 
ED and UHF results is done, one can immediately observe a 
spectacular agreement; a match of the order of 99.9% in aver- 
age has been reached for both local charge and spin densities. 
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FIG. 5: (Color online), (left column) Comparison between ED 
(small filled symbols) and UHF (large open symbols) results of the 
local charge (circles) and spin (triangles) densities as a function of 
the site index in the case of two defect plaquettes. (right column) 
Local charge and spin density snapshots, (a) 4 electrons in the 36 
site cluster (S = 1) for (U, V) = (1.5, 1). (b) 4 holes in the 64 site 
cluster (S = 2) for (U, V) = (1.5, -0.2). 



The texture of the local charge and spin densities is similar in 
each case and exhibits interesting and rich patterns, depend- 
ing on both the positions of the defect plaquettes and on the 
carrier concentration. We have considered two different situ- 
ations in case (a) and (b) of Fig. |3J in (a) the defect plaquette 
positions break the lattice (point group and translation) sym- 
metries while some remain preserved in case (b). We consid- 
ered electron and hole doped for respectively (a) and (b) and 
the same strength for U. In Fig. |5ja), one sees that the local 
spin density fluctuates very strongly from site to site and is 
found very tiny on the plaquettes sites, although V is not that 
large. In Fig[5]L)) case, one first clearly sees the appearance 
of stripes and since the GS is fully polarized, the spin and 
charge textures are identical. Note that this particular stripe 
structure is a natural consequence of the symmetry of the pla- 
quette defects on a cluster with periodic boundary conditions. 
We have checked that this peculiar stripe structure survives 
on larger equivalent clusters (same defect pattern and carrier 
density). In addition one observes a clear dominant moment 
on the plaquette sites although the strength of V is only -0.2, 
which is very small. Fig. [6] on the graphene lattice, is even 
more remarkable since a very inhomogeneous and complex 
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FIG. 6: (Color online). Local charge and spin densities in a 54-site 
hexagonal lattice as presented in FigfJ] for 4 electrons and (U, V) = 
(0.2, 0.2). The spin ground state is S = 1. 



GS is obtained even though very small values of both parame- 
ters (U = V = 0.2) were considered. Indeed, for the spin tex- 
ture, the ratio between the largest and the lowest local density 
values (around the defect plaquettes) is at least of the order 
10. More precisely, the Honeycomb lattice case is the most 
impressive. The fluctuation between the lowest and the high- 
est local densities is at least 300%. From the previous two 
figures, it is now clear that such nano-patterned or decorated 
lattices present a very rich and complex physics. Moreover, 
the UHF approach appeared to be spectacularly accurate in 
the description of the GS properties which was far beyond our 
expectations. It is important to emphasize that we have also 
checked whether non-collinear phases could be the GS by in- 
cluding explicitly the transverse field: 



(5) 



In all investigated cases, only collinear phases have been 
found. 

Let us now show that the accuracy of the UHF is not re- 
stricted to weak and intermediate coupling only. 

In order to illustrate this, we compare up to large U, both 
local charge and spin densities obtained within UHF and ED. 
For that purpose, we define the following standard deviation 



a: 



-Ed 



(6) 



where A, = £ 2 UHF /£f D and £j is either the local spin Si tZ or 
the local charge p. L . We focus on the case of Fig.Uc) (6x6 
triangular lattice with two plaquettes and filled with 4 elec- 
trons) for V = — 1 and U ranging from 1 to 12. We recall 
that for U > 1, the total spin is S = 1, thus an average 
spin per site of = In our definition, Njj — N for 
= pi whilst for the local spin we keep only the sites for 
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This restriction is introduced in or- 
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sites with tiny values of Sj lZ . In other words, negligible values 
of Si z are disregarded. However, typically, Ny is always at 
least of the order of ~ 0.8iV. The results are depicted in Fig|7] 
On clearly sees that up to U = 4, the agreement is of the order 
of 99% for both local charge and spin densities. 
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FIG. 7: (Color online). Standard deviation between UHF method 
and ED for both the local spin density and local charge as a function 
of U. The system is a 6 x 6 site cluster on the triangular lattice filled 
with 4 electrons. The plaquette potential is set to V = — 1, 

As U increases, a increases monotonously and becomes 
flatter for large U. However, even in this strong coupling 
regime (U = 12), we obtain an excellent agreement of the 
order of 95% and 90% for pi and Si >z respectively. 

One crucial and natural question which arises now is 
whether the ferromagnetism survives in the thermodynamic 
limit. For that purpose, we have computed as a function of 
the system size up to N = 84 x 84 sites the magnetic mo- 
ment per site S/N, for a given pattern on the triangular lat- 
tice. The results are depicted in Fig. [8] for various hole con- 
centrations rih = 2 — n where 1.5 < n < 1.7. Remark that 
n = 1.5 is particular since the Fermi energy coincides with 
the van Hove singularity of the host triangular lattice density 
of state (U = V = 0). We observe for that particular den- 
sity a singlet phase S = for any system size. However, the 
GS nature changes completely once the hole concentration is 
slightly decreased (n increases to ~ 1.53). Indeed, we ob- 
serve an almost finite constant S/N ~ 0.015 ratio w.r.t. 1/N. 
For example, at N = 20 x 20, a total spin GS of S = 6 
has been found. As n is further increased to 1.56, one ob- 
serves a significant jump of S/N ~ 0.03 and a similar be- 
havior with the system size, revealing a strong sensitivity to 
the hole concentration. Note that the tendency of the data to 
weakly deviate from the linear fit has two possible origins. 
First, due to commensurability, the carrier density is set in 
average, e.g. for N — 36 x 36, n = 1.52932 instead of 
1.53. Then, as the system size increases and as mentioned be- 
fore, a very small change in the hole concentration leads to 
strong S/N variations. For n ~ 1.60, one now sees a differ- 
ent finite size behavior of S/N which is a linear increase as 
a function of N. It is remarkable to notice that in the ther- 
modynamic limit, S/N coincides with that of n ~ 1.56. For 
the largest hole concentration considered n — 1.70, we find 
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FIG. 8: (Color online). Spin ground-state per site S/N as a func- 
tion of 1 /N on a nano-pattemed triangular lattice, TV being the total 
number of sites defined as N — I x I, with I up to 84. The param- 
eters are U = 1.5 and V = 1. The two vectors characterizing the 
nano-pattern are ti = (2,2) and t2 = ( — 2,2) (see Fig. [TJ. The 
electron concentration n varies from 1.50 to 1.70. Continuous lines 
correspond to linear fits of the UHF data (symbols). 



a surprising non monotonic behavior of S/N. However, in 
the thermodynamic limit, the spin GS is a singlet. In the light 
of these results, the unanticipated strong hole concentration 
sensitivity of S/N for such patterns should open interesting 
perspectives for spintronic devices. Indeed, by tuning the car- 
rier density (gate voltage, electron/hole co-doping) on such 
2D nano-patterned lattices, one can switch and manipulate the 
magnetic moment of the material. 

In conclusion, the competition/interplay between kinetic 
term, Coulomb repulsion and correlated on-site disorder leads 
to very rich and complex ferromagnetic phase diagrams on 
nano-patterned / decorated lattices. In order to illustrate the 
richness of the underlying physics, we have considered sys- 
tems of various geometries and natures. Depending on the 
carrier nature (electron/hole), its concentration, the type of lat- 
tice, the region of stability of ferromagnetism can be very ex- 
tended, and exhibits multiple phase transitions with even the 
possibility of re-entrance in some cases. The phase diagrams 
depend strongly on the chosen patterns. In addition, even in 
the case of relatively small on-sit potential, the ground state 
appears to be very inhomogeneous. Huge fluctuations of both 
the local spin and charge densities have been evidenced. We 
have also found that, in the thermodynamic limit, the ground 
state nature can change drastically even for a small variation 
of carrier density. The other significant output was to estab- 
lish that an accessible mean field type approach, namely the 
Unrestricted Hartree-Fock method, is a promising and power- 
ful tool to investigate the many body physics on large scale 
systems. An excellent but unexpected quantitative agree- 
ment, even at large U, between UHF and ED have been evi- 
denced. With the support of both theoretical based approaches 
(ab-initio/DFT) and the tremendous expertise in customizing 
nano-cluster arrays at the atomic precision ?4 , the realization of 
such nano-patterned compounds is already within our reach. 
Thus, we hope that these findings will open new routes to de- 
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sign real materials that could be good candidates for spintronic 
devices. 
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